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One-loop corrections to the Drell-Yan process in SANC 
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Abstract. Radiative corrections to the charged current Drell-Yan processes are revisited. Complete 
one-loop electroweak corrections are calculated within the automatic SANC system. Electroweak 
scheme dependence and the choice of the factorization scale are discussed. Comparisons with 
earlier calculations are presented. 
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PACS. 1 3.85.Qk Inclusive production with identified leptons, photons, or other nonhadronic 
particles; 12.15.Lk Electroweak radiative corrections 

1 Introduction the Drell-Yan process at high energy hadronic col- 
lider were calculated by several groups in the past, see 

Precision studies of the Drell-Yan process are vitally paperg rTTElEflCTMf)! and references therein. Here we 

important for high energy hadronic colliders. This pro- present the results for the corrections to the charged 

cess provides information about weak interactions and current Dre ll_ Y an process, obtained within the autom- 

contributes to the background to many of the searches atized sygtem SANC rjjrgj and some comparisons with 

for physics beyond the Standard Model. One-loop QED earlier ca i cu i at ions. Starting from the construction of 

and electroweak (EW) radiative corrections (RC) to helicity amplitudes and EW form factors, SANC per- 
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forms calculation of the process cross section and pro- As usually, we subdivide the EW RC into the vir- 
duces computer codes, which can be further used in tual (loop) ones, the ones due to soft photon emission, 
the experimental data analysis. and the ones due to hard photon emission. An auxil- 

iary parameter Q separates the soft and hard photonic 
contributions. 

2 Preliminaries and Notation 

In the automatized system [5], the virtual correc- 
Let us start with the partonic level, where we will con- tions are accessible via menu chain SANC — > EW — > 
sider interactions of free quarks (partons). The differ- Processes — ► 4 legs — > 4f — > Charged Current — > 
cntial Born- level cross section of the process fl fl' — > f P (FF). The module, loaded at the end 



d(pi) + u{p 2 ) -> l+( P4 ) + i4(p 3 ) (1) 



of this chain computes on-line the scalar form factors 
of the partonic sub- process (JIJ. The parallel module 
in the center-of- mass system of the initial quarks reads . . . fl fl' — > f f (HA) provides the relevant helic- 



j » *-,2 n/r2 ~2 ity amplitudes. For more details see Section 2.5 of the 

dJ? 4A^ C " d 2tts (s-M2,) 2 + r2,(s)M2,' SANC description 7 and the book 0. 
s=(pi+p 2 ) 2 , "=(P1-P 3 ) 2 , (2) 

1 ne real pnoton emission process 

where N c = 3 is the number of quark colors; V ud is J( Pl ) + u(pa) -> l + {Pi) + ^(Pa) + 7(ps) (3) 
the relevant element of the CKM matrix; Gf is the 



Fermi coupling constant; M w and r w are the mass 
and the width of the M^-boson, respectively. 



should be taken into account as well. Integration over 
the phase space in this case can be performed either 
(semi-)analytically or by means of a Monte Carlo in- 
tegrator. 

3 Radiative Corrections at the Partonic The nrst possibility is realized within the SANC 



Level 



environment. Now we have there two branches. The 
first one contains the complete chain of analytical in- 
In order to get a more accurate description of the pro- tegrals over the hard photon phase space. It provides 
cess we should go beyond the Born approximation and at the partonic level the double-differential distribu- 
take into account different sources of radiative correc- tion d 2 <Thard/(dc ds') and the single differential dis- 
tions. Here we will consider only EW contributions tribution d<7hard/dc, where c = cosZ(p2P4) and s' = 
to the corrections, while effects of higher order QCD (p^ + pi) 2 . The second branch provides the double- 
contributions (and mixed effects) are left beyond the differential distribution d 2 <Thard/(dc dM 2 ), where M 2 = 
scope of our study. %P3P5 which is directly related to the charged lep- 
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ton energy in the center-of-mass system of the initial 
quarks: 



- o s + mf- Ml 



(4) 



We managed also to obtain analytically the hard 
photon contribution as the single differential distribu- 
tion d<7hard/ds'. In this case we can use a system of 
reference with the z-axis along the real photon mo- 
mentum j>5 . There the integration over three angular 
variables is rather easy and we have a possibility even 
to keep all the light masses exactly. Below we give the 
expression without mass terms because of its simplic- 
ity: 



d<T h 



aid 



a 1 



1 



d,s' 



x In 



2tt s 2 s 



Is 2 + l ,2 l 



Qf 



s {s'-MiY + riMi 



1 + Qi In 



(? + $$' + 8") 



Qf 



$> {s-Mlf + rlMl 



s {s'-Mlf+TlMl 
1 [§'(§ + §')} Qi(4Q u + 5Q d } 
\Ml-s) (Ml-s')+rlMl 



(5) 



{s'-Miy + n.Mi y 

where Qi, Q u , and Qd are the charges of the charged 
lepton, up-quark and down-quark, respectively. 

The differential distributions of the tree-level ra- 
diative process d+u — > fi + +i> fl + r y were compared with 
the corresponding distributions obtained by means of 
the CompHEP package JOI- Cross section distribu- 
tions in the cosine of the outgoing muon angle and 
in the muon energy are considered. 20 bins are con- 



structed for each of the distributions. Bins in the muon 
energy are 

{n hin - 1) x 5 GeV < < n hin x 5 GeV. (6) 

The cut on the muon energy (E^ < 95 GeV) is im- 
posed in both distributions to avoid the region with 
soft photons, where CompHEP is not supposed to 
work well. The angular bins are 



(7) 



The a(Mz) electroweak scheme (realized according to 
the CompHEP conventions) was used. An agreement 
was found as can be seen from Tabled 

For the two choices of variables we have simple an- 
alytical expressions for the corresponding soft photon 
contributions. The infrared singularities in them are 
regularized by the auxiliary photon mass. The energy 
of a soft photon is limited from above by a cut in the 
integral either over s' or over M%. 

In order to have the possibility to impose exper- 
imental cuts and event selection procedures of any 
kind, we can use a Monte Carlo integration routine 
based on the Vegas algorithm 11 . In this case we 
perform a 4(6)-fold numerical integration to get the 
hard photon contribution to the partonic (hadronic) 
cross section. To get the total EW correction we add 
also the contributions of the soft photon emission and 
the ones of the virtual loops. The cancellation of the 
dependence on the auxiliary parameter ui in the sum 
is observed numerically. 
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Using the splitting of the VF-boson propagators in theorem j!3U14| . But if the experiment is measuring 

the case of real photon emission off the virtual W, we the energy of the charged lepton without summing it 

separate the contributions of the initial state radia- with the energies of accompanying collinear photons 

tion, the final state one, and their interference in a (typical for muons), the logarithms with the lepton 

gauge invariant way |12j . The splitting is introduced mass singularity remain in the result and give a consid- 

by the following formula: erable numerical contribution. Re-summation of these 

i ^ logs in higher orders was discussed in Refs. |ToTll6lll7| . 



s - (m w - ir w y §' - {m w - ir w y 
if i 



{s-s')\s'-{M w -ir w ) 2 

i \ 3.1 Treatment of Quark Mass Singularities 

(8) 



s-{M w -ir w f i 

In the center-of-mass system (s - §') = 2pgVI. The Qne-loop radiative corrections contain terms propor- 

fixed W-width scheme is used here and in what fol- tional tQ the logarithms of the quark masseSj \ n (S/m 2 u d ). 

l° ws - They come from the initial state radiation contribu- 

In the course of calculations of the O {a) correc- tions inc i uding hard, soft and virtual photon emission, 

tions we met the so-called on-shell singularities, which Such Mtial state mass sin g U i ari ties are well known, 

appear in the form of ln(J - M 2 , + ie). As was shown for instance; in the process of e+e - annihilation. But 

in detail in Ref. 0, they can be regularized by the in the case of hadron collisions these logs have been 

VT-width. already effectively taken into account in the parton 



ln(s' - Ml + ie) -» ln(s' - M? v + iM w r w ). (9) 



density functions (PDF's). In fact, in the procedure of 



PDF's extraction from the experimental data, QED 

In the analytical formulae for radiative corrections rad i ative corrections to the quark line have not been 

one can find logarithms with quark and lepton mass systematically subtracted. Therefore the present PDF's 

singularities: effectively include not only the QCD evolution but 

s s s also the QED one. Moreover, it is known that the lead- 

in — 5-, In — 5-, In — 7T . (10) 
mf rat. mi 

1 u a ing log behaviors of the QED and QCD DGLAP evo- 

In the experimental set-up with calorimetric registra- lution of quark density functions arc similar (propor- 
tion of the final state charged particles (typical for tkmal to each othcr) ^ go Qne getg evolution of FBPs 
electrons), the lepton mass singularity cancels out in with an effective coup i ing constant 

the result for the correction to an observable cross see- 
off Qi 

tion in accordance with the Kinoshita-Lee-Nauenberg a s ~ a s + ~Q~^ a ' 
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where a s is the strong coupling constant, a is the fine where Qi and rrii denote the charge and the mass of 

structure constant, Qi is the quark charge, and Cf the given quark; M is the factorization scale; (To(^) 

is the QCD color factor. The nontrivial difference be- is the cross section at the partonic level with the re- 

tween the QED evolution and the QCD one starts to duced value of the quark momentum: pi —* £,iPi- The 

appear in higher orders, and the corresponding numer- subtracted partonic cross section with O (a) correc- 

ical effect is small compared to the remaining QCD un- tions is given by 

certainties in PDF's 08 , 19 20 ,21 . The best approach 

of s = &! - S us . (13) 

to the whole problem would be to re-analyze all the 

experimental DIS data taking into account QED cor- Then it can be convoluted with PDF's as shown below 

rections to the quark line at least at the next -to- in Eq. (|16(l . 

leading order. But for the present moment we can limit But there is an alternative way to perform the sub- 

ourselves with an application of a certain subtraction traction. Really, to avoid the double counting of the 

scheme to the QED part of the radiative corrections quark mass singularities, we can leave them in the cor- 

for the process under consideration. We will use here rected cross section, but remove from the PDF's: 



the MS scheme [22], the DIS scheme can be used as r 1 a / 



well. This allows to avoid the double counting of the 
initial quark mass singularities contained in our result 



q(x, M 2 ) — q(x, M 2 ) — ^ q (l,M*\±Q 



1 — z I V m 2 



ln( — ) -21n(l-js)-l[ 

r2\ 



for the corrections to the free quark cross section and = <l{x, M ) — Aq, (14) 

the ones contained in the corresponding PDF. The lat- 

where q(x, M 2 ) can be taken directly from the exist- 

ter should be also taken in the same scheme with the 

ing PDF's in the MS scheme (see Ref. [3] for the corre- 

samc factorization scale. 

sponding formula in the DIS scheme). It can be shown 

In fact, using the initial condition for the non- 

analytically (see i.e. Ref. PJ), that this procedure is 

singlet NLO QED quark structure function, which co- 
equivalent to the subtraction from the cross section, 

incides with the QCD one with the trivial substitution 

and that it really removes (hides) the dependence on 

C^tts — > Q 2 a, see Ref. (22], one gets the following ex- 

the quark masses. The advantage of the last approach 

pression for the terms to be subtracted from the full 

is that it can be used regardless of the way to repre- 

calculation with massive quarks: 

sent the partonic cross section: it can be kept even in 

l 



i=l,2 



21n(l-&)-l 



M 2 the completely differential form. 



In 

mf 



<5o(&)i (12) 1/2 , r . 



+ 



The natural choices of the factorization scale are 
M 2 = (when the returning to the VF-resonance 
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is allowed by kinematic cuts) and M 2 = s = x\x 2 s. where the parton densities with bars mean the ones 
Variations with respect to the choice should be stud- modified by the subtraction of the quark mass singu- 



ied. 



larities; the step function <9(c, E^) defines the phase 



In order to avoid the appearance of spurious higher space domain corresponding to the given event selec- 
order terms for the case of subtraction from PDF's, we tion procedure. The partonic cross section is taken 
suggest to apply a procedure of linearization. Schemat- in the center-of-mass reference frame of the initial 



ically it can be represented as follows: 

qi(x!,M 2 ) x q 2 {x 2 ,M 2 ) xoi = [ qi ( Xl , M 2 ) - Aq x ] 
x [q 2 (x 2 , M 2 ) - Aq 2 ] x ((TBorn + a a ) 

-> qi(xi,M 2 ) X q 2 (x 2 ,M 2 ) X (TBorn 

+q 1 (x 1 ,M 2 ) x q 2 (x 2 ,M 2 ) x a a (15) 
-[ qi (x u M 2 ) X Aq 2 +q 2 (x 2 ,M 2 ) x A qi ] X 

CBorn, 

where o^om and a a denote the Born-level partonic 
cross section and the O (a) RC contribution to it, re- 
spectively. Without the linearization procedure, terms 
with quark mass singularities would remain in the 
O (a 2 ) contribution to the cross section. 



quarks, where the cosine of the muon scattering angle, 
c, and the muon energy, E^, are defined. The transfor- 
mation into the observable variables c and involves 



the Jacobian: 



„ dc dE u 

J = TT 



ix x x 2 



a 2 (l + c) 



dc dE u 



a = X\ + x 2 — c(xi — x 2 ), 



Xi[a + x 2 (l + cj] 

- 1 \ 2xi 
c = 1 - (1 - c) — , 



SXiX 2 , 



Efi — E M 



E/, — — , 
2 



(17) 



1-c 2 ' 

An analogous formula can be written for any other 
choice of a differential distribution as well as for the 
total cross section. 



4 Radiative Corrections to Hadronic 
processes 

The double-differential cross section of the Drell-Yan 
process can be obtained from the convolution of the 
partonic cross section with the quark density func- 
tions: 

i+vX, 



dc dE n 



/ / dxi dx 2 q 1 (x 1 ,M 2 ) 



xq 2 (x 2 ,M 2 ) 

dc ah,, 



J9{c,EJ, (16) 



5 Numerical Results and Conclusions 

For numerical evaluations we take the same set of 
input parameters as the one given by Eq. (4.1) of 
Ref. U2|. In Table|21we present the results for the total 
cross section 1 of the process u + d — > vi + Z + (+7). For 
the Born-level cross section we completely (in all listed 
digits) agree with the numbers given in Ref. 0. The 
third line shows radiative corrections in percent be- 
fore the subtraction of quark mass singularities. These 
1 factor |Vud| 2 has been dropped in the sake of compar- 
ison with Ref. 0. 
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numbers were received directly from the SANC system case we are going to perform further studies of this 

for the Gf EW scheme. Starting from the fourth line effect. 

we use the treatment of the EW scheme 2 , which has m 1 1 m ^ j j r ii 1 j 

1 a ble|4| represents the dependence of the hadromc 

been adopted Ref. 151. The results for the radiative „ v ,. , r , , , 

r 1-1 Drell-Yan cross section on the values of the quark 

corrections with MS subtraction (with factorizations ., , , . ,, , 

v masses with and without the subtraction procedure. 

scale being equal to Mw) are also in a fair agree- T , c u 4-1 ± c 

^ *" ' a the conditions are as follows: the center-of-mass en- 

ment. The small deviations there can be due to de- • nn r~\ u n , ^ - r 

ergy is 200 GeV, all events with the invariant mass of 

tails in the treatment of EW scheme with respect in- , , . , , , , , rrx „ , T 

v the neutrino and charged lepton pair above y40 GeV 

duced higher order effects. Huge positive corrections , , ji.xi.Dii 

are accepted, cto denotes the Born-level cross section 

in the case without subtraction of quark mass singu- , , . , . nrr ^njT ± r nnn rrn 

H to obtained using the CTEQ4L set of PDF s |23j. <ti, 

larities above the VK-peak are due to the initial state MS(o-) ms(o) , ms(?)/,. s , , r 

^ fi]^ , (Jj^ , and cr, w (lin.J stand for the cross 

radiation which provides the radiative return to the . in , „,„ , , , ml , , , 

1 sections with one-loop EW RC included, ihe double 

VK-resonance. c ^ 1 1 • 

counting of the quark mass singularities 111 o\ is not 

The effect of EW scheme dependence is illustrated removed. The MS procedure |T3J) is applied to the par- 
by Table El Results for the total partonic cross sec- tonic cross section in the computation of crf S(<T) . Val- 
tion at the Born and O {a) levels are given for two ues of crf S(?) and af S(9) (lin.) are computed by con- 
EW schemes. At the Born level the 7.3% difference volution of the quark (parton) density function mod- 
appears just due to the difference in the definition of ified according to Eq. (0! with the full (including 
EW constants in the G F and in the a(0) schemes. As quark mass singularities) partonic cross section. The 
it should be the difference between the corrected cross linearization procedure JT5J was adopted for erf S{q) (lin. 
sections is less than the one at the Born level. But still in addition. One can see that the numerical effect of 
it is large and comparable with the ordered precision linearization for the given set-up is small (but visible), 
of the calculation. Certainly, usage of the a(0) is not The two approaches to remove the double counting 
well motivated for the given energy range. And the S ive vei T dose results as it should be. 
difference Si gives only an upper estimate of the un- For an internal test of our calculations, a corn- 
certainty due to the EW scheme dependence. In any par ison of the results produced by our Monte Carlo 

(MC) and semi-analytical (SA) codes for the descrip- 



2 In G' F scheme we assigned the following one-loop value tion of hard photon contributions was performed. The 
of the coupling constant standing at the photon vertices: results are presented in Tablc|5J where the correspond- 
olqed ~ 1/132.544. ing contributions to the proton-proton cross section 
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at 14 TeV center-of-mass energy are given. The condi- agreement with the ones published earlier in Ref. 
tions and the input parameters were the taken as the The corresponding computer codes in analytical (FORM) 



ones used in Ref. |25| : 



G p 


= 1.16637 x 10- 5 GeV" 2 : 






a(0) 


= 1/137.03599911, 


a s 


= 0.1187, 


M w 


= 80.425 GeV, 


r 

1 w 


= 2.124 GeV, 


M z 


= 91.1867 GeV, 


r z 


= 2.4952 GeV, 


M H 


= 150 GeV, 


m t 


= 174.17 GeV, 


m u 


= nid = 66 MeV, 


m c 


= 1.55 GeV, 


m s 


= 150 MeV, 


m b 


= 4.5 GeV, 


\v ud \ 


= \V CS \ = 0.975, 


\v us \ 


= \V cd \= 0.222 



and numerical (FORTRAN) formats are available from 
SANC They can be used as a part of a more 
general computer program (like a Monte Carlo event 
generator) to describe the Drell-Yan process in real- 
istic conditions. Further comparison at the hadronic 
level with analogous calculations of other groups is in 
progress |2*5] . 
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6 Tables 





bins 


in E,, 


bins in c 


bin 


SANC 


CompHEP 


SANC 


CompHEP 


1 


0.0006(1) 


0.0006(1) 


0.5867(1) 


0.5869(2) 


2 


0.0010(1) 


0.0010(1) 


0.3538(1) 


0.3537(1) 


3 


0.0023(1) 


0.0023(1) 


0.1857(1) 


0.1858(1) 


4 


0.0452(1) 


0.0451(1) 


0.1111(1) 


0.1111(1) 


5 


0.0569(1) 


0.0569(1) 


0.0740(1) 


0.0741(1) 


6 


0.0549(1) 


0.0549(1) 


0.0541(1) 


0.0541(1) 


7 


0.0546(1) 


0.0546(1) 


0.0427(1) 


0.0428(1) 


8 


0.0563(1) 


0.0563(1) 


0.0360(1) 


0.0361(1) 


9 


0.0603(1) 


0.0603(1) 


0.0321(1) 


0.0321(1) 


10 


0.0667(1) 


0.0666(1) 


0.0298(1) 


0.0297(1) 


11 


0.0755(1) 


0.0755(1) 


0.0287(1) 


0.0287(1) 


12 


0.0868(1) 


0.0867(1) 


0.0284(1) 


0.0285(1) 


13 


0.1008(1) 


0.1008(1) 


0.0286(1) 


0.0287(1) 


14 


0.1175(1) 


0.1175(1) 


0.0292(1) 


0.0292(1) 


15 


0.1372(1) 


0.1372(1) 


0.0299(1) 


0.0299(1) 


16 


0.1605(1) 


0.1605(1) 


0.0302(1) 


0.0301(1) 


17 


0.1881(1) 


0.1881(1) 


0.0294(1) 


0.0293(1) 


18 


0.2227(1) 


0.2227(1) 


0.0263(1) 


0.0263(1) 


19 


0.2739(1) 


0.2740(2) 


0.0187(1) 


0.0187(1) 


20 


0.0 


0.0 


0.0059(1) 


0.0059(1) 



Table 1. Bin by bin comparison of differential distributions for the process d + u — > /i + + v^, + 7. 
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VI/GeV 


40 


80 


120 


200 


500 


1000 


2000 


o-o/pb 


2.646 


7991.4 


8.906 


1.388 


0.165 


0.0396 


0.00979 


S/%, full, Gp 


-1.70 


-7.62 


89.9 


125.3 


155.9 


166.9 


173.0 


S/%, full, G' F 


-1.76 


-7.87 


92.9 


129.6 


161.2 


172.5 


178.9 


S/%, MS(s), G' F 


0.56 


2.48 


-17.2 


-16.0 


-17.9 


-23.2 


-31.5 


S/%, MS(M W ), G' F 


0.73 


2.48 


-12.9 


-3.2 


11.6 


18.8 


22.8 


<V%, 


0.7 


2.42 


-12.9 


-3.3 


12 


19 


23 



Table 2. The total lowest-order partonic cross section <ro in the Gf EW scheme and the corresponding relative one-loop 
correction 5. 



v/I/GeV 


40 


80 


120 


200 


500 


1000 


2000 


ao/pb, [Gf] 


2.646 


7991.4 


8.906 


1.388 


0.165 


0.0396 


0.00979 


ao/pb, [o(0)] 


2.454 


7410.2 


8.258 


1.287 


0.153 


0.0368 


0.00908 


So/% (diff) 


-7.3 


-7.3 


-7.3 


-7.3 


-7.3 


-7.3 


-7.3 


ffi/pb, MS (M ff ), [G F ] 


2.665 


8183.2 


7.796 


1.345 


0.183 


0.0467 


0.01195 


ai/pb, MS (M w ), [a(0)] 


2.617 


8029.5 


7.721 


1.324 


0.179 


0.0455 


0.01162 


Si/% (diff) 


-1.8 


-2.0 


-0.5 


-1.5 


-2.6 


-3.1 


-3.3 



Table 3. The total partonic cross section in the Gf and a(0) EW schemes. 





(jo [pb] 


(71 


MS(o-) 
CT 1 


MSfg) 


MSfoln- \ 
a l w (litl.) 


m u =m d = 4.85 MeV 


2.5577(1) 


2.4795(3) 


2.5724(3) 


2.5704(3) 


2.5729(3) 


m u =m d = 48.5 MeV 


2.5577(1) 


2.4992(3) 


2.5724(3) 


2.5713(3) 


2.5727(3) 


m u =m d = 485. MeV 


2.5577(1) 


2.5190(3) 


2.5724(3) 


2.5719(3) 


2.5726(3) 



Table 4. The tree level and corrected hadronic Drell-Yan cross section different values of the light quark masses. 
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p T [GeV] 


> 25 


> 50 


> 100 


> 200 


> 500 


> 1000 


CTBorn [pb] 


2112.22(2) 


13.1507(2) 


0.94506(1) 


0.115106(1) 


0.00548132(6) 


0.000262108(3) 




cD = 0.01 GeV 


fohard [%] 


27.52(2) 
27.54(2) 


34.02(2) 
34.02(2) 


43.88(2) 
43.87(2) 


51.22(2) 
51.21(2) 


59.67(2) 
59.68(2) 


65.22(2) 
65.22(2) 


5<T%f [%} 

<5<rfi [%] 


-2.69(2) 
-2.68(2) 


-3.83(2) 
-3.83(2) 


-6.67(2) 
-6.68(2) 


-11.77(2) 
-11.78(2) 


-22.28(2) 
-22.27(2) 


-33.36(2) 
-33.36(2) 




Q = 0.001 GeV 


[%] 
[%] 


36.85(2) 
36.88(2) 


44.35(2) 
44.37(2) 


55.70(2) 
55.70(2) 


64.25(2) 
64.26(2) 


74.14(2) 
74.15(2) 


80.72(2) 
80.72(2) 


5^ t ° [%] 
<5afi [%] 


-2.70(2) 
-2.67(2) 


-3.85(2) 
-3.83(2) 


-6.69(2) 
-6.69(2) 


-11.76(2) 
-11.76(2) 


-22.29(2) 
-22.28(2) 


-33.35(2) 
-33.35(2) 



Table 5. The hadronic Drell-Yan cross section in the Born approximation and the hard photon contributions to it for 
different values of the cut on the muon transverse momentum. 



